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We demonstrate that the zero-temperature conductance of the Anderson model can be calculated 
within the Landauer formalism combined with static density functional theory (DFT). The pro- 
posed approximate functional is based on finite-temperature DFT and yields the exact Kohn-Sham 
potential at the particle-hole symmetric point. Furthermore, in the limit of zero temperature it 
correctly exhibits a derivative discontinuity which is shown to be essential to reproduce the conduc- 
tance plateau. On the other hand, at the Kondo temperature the exact Kohn-Sham conductance 
overestimates the real one by an order of magnitude. To understand the failure of DFT we resort 
to its time-dependent version and conclude that the suppression of the Kondo resonance must be 
attibuted to dynamical exchange-correlation corrections. 
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Despite the many successes of density functional the- 
ory (DFT) [l[ in the description of the electronic struc- 
ture of many- electron systems, the treatment of strongly 
correlated systems is a notoriously difficult challenge. 
However, the fundamental theorems of static and also 
time-dependent (TD) DFT are exact. Therefore in 
principle these systems should be accessible and the fail- 
ure to describe them within (TD)DFT has to be ac- 
credited to shortcomings of the approximations for the 
exchange-correlation (xc) functional used in practice. 

An exact property of the static xc functional is the 
derivative discontinuity at integer occupancy Re- 
cently the crucial importance of this property in the de- 
scription of strong correlations has been recognized in, 
e.g., the Hubbard model 0, [f|, gaps in strongly corre- 
lated molecules 6|, Coulomb blockade (CB) [7[ and the 
Mott transition 8[. 

For non-equilibrium situations, the development of TD 
xc functionals is still at its infancy. Typically one uses 
the so-called adiabatic approximation by importing func- 
tionals constructed for static DFT to the time domain. 
In the context of quantum transport, one of the most 
debated questions is whether (or under which circum- 
stances) static DFT within the Non-Equilibrium Green's 
Function (NEGF) approach can yield accurate conduc- 
tances. In a seminal paper, Schmitteckert and Evers @ 
provided strong numerical evidence that the Kohn-Sham 
(KS) conductance of single channel correlated model 
junctions is in good agreement with the exact one. Later, 
the accuracy of the KS conductance has been explained 
and assessed using the Friedel sum rule [l(j ■ In general, 
however, the DFT+NEGF formalism is incomplete. A 
proper TDDFT treatment [ll|, [l2| leads to dynamical xc 



corrections [13( , the quantitative importance of which has 
yet to be clarified. 

In this Letter we show that the derivative discontinuity 
is the necessary and sufficient ingredient to describe the 
Kondo effect [14j in the zero temperature conductance 
using DFT+NEGF. We propose an analytic KS poten- 
tial for the Anderson model based on finite-temperature 
DFT for an isolated impurity. For this approximate 
potential, the derivative discontinuity and, as a conse- 
quence also the conductance plateau, emerge naturally 
in the zero-temperature limit. At finite temperatures, 
instead, we demonstrate that the DFT+NEGF approach 
is not sufficient. Although our potential is exact at the 
particle-hole (ph) symmetric point, at the Kondo tem- 
perature the KS and the exact conductances differ by 
almost one order of magnitude. We thus give a first ex- 
plicit example where the exact dynamical xc corrections 
can be quantified and show that they can be as large as 
the quantum of conductance. 

We consider the Anderson model for a single non- 
magnetic impurity attached to two non-interacting leads. 
The Hamiltonian of this system is given by 



H = Hi+ J2 



H 



T ■ 



(1) 



a=L,R 



Here, H a = - ( v 4+ia,<xCicx,a + h.c.) describes, 

in standard notation, the tight-binding lead a = L, R 
while H T = -J2 a ,a( V ^ c[ aa d a + h.c.) accounts for 
the (symmetric) coupling between the impurity and the 
leads. In this work we focus on the half-filled system and 
take Mink C V. In this parameter range the only relevant 
energy scale for electron tunneling is 7 = 2V l f nk /V (wide- 



band limit). The impurity Hamiltonian is 

Hj = Voh + Uh^n± 



(2) 



where vq is the on-site energy (or gate voltage), U is the 
charging energy, n a = d\d a is the number operator for 
electrons of spin a at the impurity and n = + nj_. 

For a DFT treatment of the problem, the first task is 
to construct an approximation to the xc potential. Here 
we propose a KS potential based on an isolated impurity 
with Hamiltonian Hj in contact with a thermal bath at 
inverse temperature j3 and chemical potential /z [l6| . The 
non-interacting KS Hamiltonian then reads 



m 



v s n 



(3) 



For both Hamiltonians Hi and Hj, the eigenstates for 
electron occupation zero, one, and two are, respectively, 
|0), | t)i I I), an d | ti) with eigenvalues 0, v , v , and 
2v + U for Hi and 0, v s , v s , and 2v s for Hj. The density 
of the interacting impurity is 



n = Tr < e 



i)/Z 



(4) 



where Z = Tr |e ^™)| is the grand-canonical par- 

tition function. Equation (j4|) depends only on vq = vq — u 
and the function n(vo) can be inverted explicitly as 



Mn) = - U -? ln { —n J ' 

(5) 

with Sn = n—1. The Hartree-exchange-corrclation (Hxc) 
part of the KS potential can then be obtained from 

vnxc{n) = v s (n) ~ v (n) = — + g{n - 1) (6) 



where g(x) = ^ + ^ In I — * 1+x I . This is 

an odd function of its argument, g(— x) = —g(x), and 
therefore vhxc(« = 1) = ¥ f° r a ll temperatures. We 
note in passing that Eq. © can also be obtained by 
differentiation of the Hxc part of the grand canonical 
potential with respect to n. The Hxc potential ^ is 
shown in the left panel of Fig. [T] for different values of 
the temperature T = 1/(3. In the limit T — > it becomes 
a simple step function with a step of height U at n = 1. 
Thus, the T = discontinuity of the xc potential emerges 
naturally from our grand canonical DFT treatment. 

For the Anderson Hamiltonian (JXJ) we make a local 
approximation and assume that the KS potential van- 
ishes in the leads while on the impurity is given by 
v s (n) = vq + vn xc (n). We expect this approximation 
to be accurate for U/j 3> 1. The KS Hamiltonian then 
reads H s = Hf + H L + H R + H T with Hj from Eq. ©. 
Using standard embedding techniques the solution of the 
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FIG. 1: Left panel: Hxc potential of Eq. ([6]) for different 
temperatures T = 1//3. Right panel: Hartree and Hxc self- 
consistent impurity density in comparison to exact results for 
T = and 7 = 10~ 2 . Energies are given in units of U. 



KS equations are reduced to the self-consistent solution 
of a single equation for the density n = no 



n = - 

IT 



dw 



7 



(w - v s (no)) 2 +r 



r//sM (7) 



where fp(u>) = 1/(1 + e 13 ^^^) is the Fermi distribu- 
tion. For low temperatures, the Hxc potential varies very 
rapidly for no ~ 1, i.e., in the regime in which we are in- 
terested (see below) . Therefore, instead of solving for no 
it is numerically much more advantageous to express no 
in the l.h.s. of Eq. ([7]) in terms of v s and solve for v s . 

The resulting self-consistent solution of the density as 
a function of vq exhibits a plateau between —U and 
with value no = 1 (sec right panel of Fig. QJ. This is in 
agreement with the exact solution of the Anderson model 
17| . Using the DFT+NEGF formalism we then proceed 
to calculate the conductance G from 



where 



G_ 
Go 



dui 



(u - v s (n )) 2 + 7 2 



(8) 



(9) 



is the zero-bias KS transmission function and Go = 1/n 
is the quantum of conductance. 

In the left panel of Fig. [2] we plot G as a function of i>o 
using the Hxc potential of Eq. © for T = ICU 2 Tk, with 
the Kondo temperature T K = VW-/ exp ( ™°( v °+ u ) 



4 7 u 



18| . For comparison we also report the results within 
the Hartree approximation, i.e., for v s = vo + Uno/2, and 
the accurate Numerical Renormalization Group (NRG) 
results from Ref. [3 The first remarkable feature is 
the plateau in the Hxc conductance with value Go in the 




FIG. 2: Hartree, Hxc and NRG [H conductances versus i; FIG ' 3: NRG conductances from Ref. [stars(blue)] and 

for two different temperatures. The parameters are: M = 0, Ref ' [solid(black)] against the exact Hxc conductance 

7 = 7T x I(T 2 [T K = 6.6 x 1(T 7 ]. All energies are given in [dashed(red)] at the ph symmetric point versus temperature, 
units of U. 



range — U < Vq < 0, in good agreement with the exact 
result. This is a direct consequence of the discontinuity 
of the proposed vrxc at n = 1 which guarantees that v s 
is pinned to the Fermi energy. In contrast, any continu- 
ous potential, like the Hartree potential, cannot capture 
the plateau. Physically, the plateau cannot be explained 
solely in terms of CB, according to which G should be 
peaked at the end-points of the plateau and should in- 
stead be small at the ph symmetric point vq = —U/2. 
The value G = Go at Vq = — U/2 is due to the formation 
of a resonant many-body singlet state between the spin of 
the impurity electron and the spin of the screening cloud, 
i.e., the Kondo effect. The first important finding of this 
Letter is therefore that the Kondo effect in G(T = 0) 
is within reach of the DFT+NEGF approach provided 
that the approximate i>hxc has the proper discontinuity 
at integer n. 

The theoretical possibility of describing the conduc- 
tance in the Kondo regime within DFT+NEGF can be 
understood in two ways. (1) For proportional coupling 
the Meir-Wingreen formula [20j at T = yields 

^=7 2 |g(M)| 2 7 ' Im ™ ] (10) 
Go 7 

where = (u> — vq — + 17) is the impurity 

Green's function and £ is the many-body self-energy. 
Since quasi-particles at the Fermi energy have an in- 
finitely long life-time then Im[£(//)] = 0. From Eq. (fT0|) 
we thus see that it is possible to reproduce the exact 
conductance in a KS system with v s = vq + Re[S(/ti)]. 
(2) From the Friedel sum rule we know that at T = 
the conductance of the Anderson model is completely 
determined by the density at the impurity [lOj, i.e., 
G = G(no). Since exact DFT yields the exact density 
then it must also yield the exact conductance. Note that 



the two explanations above are complementary but not 
equivalent since the exact KS potential is not rigorously 
zero in the neighborhood of the impurity. It is also worth 
stressing that the equality between the exact and KS con- 
ductances does not follow from an equality between the 
corresponding spectral functions. The latter are com- 
pletely different, although the Kondo peak (of width Tk) 
and the KS peak (of width 7) both occur at the Fermi 
energy. 

At finite temperature G does not depend on no alone. 
In the right panel of Fig. [5] we display the NRG, Hxc 
and Hartree conductances as a function of vq for T = 
6 • 10 4 Tk- The failure of the Hxc and Hartree approxi- 
mation in reproducing both the dip at vq = —U/2 as well 
as the CB side peaks is evident. In Fig. [3]we compare the 
Hxc conductance with the exact conductance [19j, [21[ for 
different temperatures at the ph symmetric point. While 
the Kondo peak is strongly suppressed for T > 10 Tk, the 
DFT+NEGF approach predicts a conductance Go up to 
much higher temperatures. For T > U , thermal fluctua- 
tions destroy the CB and the KS and exact conductances 
merge and approach zero for T — > 00. 

Is it possible to reproduce the temperature behavior of 
the exact G by improving the quality of the KS potential? 
The answer to this question is negative since at the ph 
symmetric point our Hxc potential is exact. For vq = 
— U/2 the Hamiltonian in Eq. ([1} is invariant under the 
transformation c^q, — > (— l) l+1 c| Qcr and d a — > a\, and 
hence the density is unity for all sites at all temperatures. 
The KS potential which reproduces this density is zero 
everywhere (in leads and impurity) for all T and Mink- 

If one cannot reproduce the exact G by using the exact 
KS potential we must conclude that the DFT+NEGF 
approach is not sufficient. To trace back the origin of the 
problem we resort to TDDFT [22[ and its lattice version 
[231 ] . In TDDFT it is possible to reproduce the exact 



4 



TD longitudinal current in a system of non-interacting 
electrons. In Ref. [ll| it was shown that if a steady-state 
is reached in the long-time limit then the current is given 
in terms of a Landauer-like formula with KS bias V a + 
V a ,xci where V a is the physical bias and V a ^ c is the xc 
correction. Then, to first order the current becomes (see 
also Ref. Il3[ ) 



I = Go{Vl-V r +Vl, xc -V r ,: 



(11) 



Consequently, the finite temperature conductance G = 
I/{Vr - V L ) coincides with that of the DFT+NEGF ap- 
proach in Eq. ((5J) only provided that V a . xc = 0. The 
exact expression for the xc bias is (to first order in V a ) 



V a 



lim 



f xc {ia,r)Sn r 



(12) 



where the sum runs over all sites r of the model, Sn r is 
the first-order density change in site r and / xc (r, r') is 
the zero-frequency xc kernel of TDDFT. In most com- 
monly used local approximations / xc (r,r') oc 8 rr > and 
hence Vq jXC vanishes since 5n r = deep inside the leads; 
as a result, it is often neglected altogether. The relevance 
of the dynamical xc correction has long been debated 
0, [24j]. The second important finding of this Letter is 
that at finite temperature the dynamical xc corrections 
are absolutely essential. They are not only important for 
the correct suppression of the Kondo peak but also for a 
quantitivc description of the CB side peaks. 

For the Anderson model we can find an explicit form 
of the dynamical xc correction to G in terms of the xc 
kernel. The linear density change at site r is 



8n r = Prr'VL 



r'£L 



^ ^ Prr' 
r'eR 



V R + P rQ ^6n Q (13) 



where P is the zero-frequency polarization. In linear re- 
sponse TDDFT P can be calculated from the xc kernel as 
P = Pq + Pq/xg-Pj Pq being the non-interacting polariza- 
tion. Substituting Eq. (|T3)) into Eq. (|12p and exploiting 
the symmetry of the L and R leads we find from Eq. (fTTj) 
that the exact conductance at any temperature reads 



| = -(i + o sc )|^r(4 

Here T(ui) is the KS transmission of Eq. © and 



Q xc = lim V V(-l) e 

l— »OG L ' L * 

r 



) fx C (ia,r)P rr ' 



(14) 



(15) 



with e a (r') = if r' € a and 1 otherwise. At T = we 
have Q xc = but for T > T K it must be Q xc > -1. The 
quantity 1 + Q xc can be interpreted as a KS dielectric 
function as it measures the ratio between the KS and the 
Hartree screening. 



In conclusion we proposed a finite-temperature DFT 
scheme to construct approximate xc functionals for cor- 
related systems weakly connected to leads. The resulting 
KS potential in the zero-temperature limit exhibits a dis- 
continuity at integer number of particles, as it should. We 
showed that the discontinuity is essential to reproduce 
the conductance plateau due to the Kondo effect within 
the DFT+NEGF approach [25| . For temperatures larger 
than 7k, however, the exact KS conductance overesti- 
mates the exact conductance by an order of magnitude. 
We traced back the origin of this problem to the lack of 
dynamical xc corrections which we expect to be relevant 
not only in this context but also in the description of 
finite temperature and finite bias transport experiments 
of weakly coupled molecular junctions. 
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